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1. Introduction 

In this work we analyse two sequential urn models with two types of balls. At the beginning, the urn contains n 
white and m black balls, with n,m > 1. Associated to the urn model are two sequences of positive real numbers 
A = (a„)„ e N and B = (/? m ) m gN- At every step, we draw a ball from the urn according to the number of white 
and black balls present in the urn with respect to the sequences (A, B), subject to the two models defined below. The 
choosen ball is discarded and the sampling procedure continues until one type of balls is completely drawn. 



ABSTRACT. In this work we discuss two urn models with general weight sequences (A, B) associated to them, A = 
(a„)„ 6 M and B = (/3 m ) m <=N, generalizing two well known Polya-Eggenberger urn models, namely the so-called sampling 
without replacement urn model and the OK Corral urn model. We derive simple explicit expressions for the distribution 
of the number of white balls, when all black have been drawn, and obtain as a byproduct the corresponding results for the 
Polya-Eggenberger urn models. Moreover, we show that the sampling without replacement urn models and the OK Corral urn 
models with general weights are dual to each other in a certain sense. We discuss extensions to higher dimensional sampling 
— without replacement and OK Corral urn models, respectively, where we also obtain explicit results for the probability mass 
functions, and also an analog of the dualitiy relation. Finally, we derive limit laws for a special choice of the weight sequences. 

o 

S 

O 

Urn model I (Sampling with replacement with general weights). Assume that n white and m black balls are contained 
in the urn, with arbitrary n, to e N. A white ball is drawn with probability a n / (a n + f3 m ), and a black ball is drawn 
with probability f3 m / (a n + (3 m ). Additionally, we assume for urn model I that an. = /?o = 0. 

Urn model II (OK Corral urn model with general weights). For arbitrary n, m e N assume that n white and to black 
balls are contained in the urn. A white ball is drawn with probability f) m / (a n + f} m ), and a black ball is drawn with 
probability a n / (a n + /3 m ) . 

The absorbing states, i.e. the points where the evolution of the urn models stop, are given for both urn models by the 
positive lattice points on the the coordinate axes {(0, n) \ n > l}u{(m, 0) | n > 1}. These two urn models generalize 
two famous Polya-Eggenberger urn models with two types of balls, namely the classical sampling without replacement 
(I), and the so-called OK-Corral urn model (II), described in detail below. We are interested in the distribution of the 
random variable X m , n counting the number of white balls, when all black balls have been drawn. For both urn models 
we will derive explicit expression for the probability mass function ¥{X n ^ m — k} of the random variable X n , m . 
Moreover, we will describe a certain duality between the sampling without replacement urn model and the OK Corral 
urn model using an appropriate transformation of weight sequences. 

> 

1.1. Weighted lattice paths and urn models. Combinatorially, we can describe the evolution of the urn models I 
and II by weighted lattice paths. If the urn contains n white balls and to black balls and we select a white ball, 
then this corresponds to a step from (to, n) to (m, n — 1), to which the weight a n / (a n + (3 m ) is associated for urn 
model I, and (3 m /(a n + (3 m ) for urn model II. Converserly, if we select a black ball, this corresponds to a step from 
(m, n) to (m — 1, n), to which the weight (3 m / (a„ + (3 m ) is associated for urn model I, and a n /(a n + f3 m ) for urn 
model II. The weight of a path after t successive draws consists of the product of the weight of every step. By this 
correspondence, the probability of starting at (m, n) and ending at (j, k) is equal to the sum of the weights of all 
possible paths starting at state (to, n) and ending at state (j, k), for j, k e N . The expressions for the probability, 
obtained by this combinatorial model are, although exact, very involved and not useful for i.e. the calculation of the 
moments of X, hm , or for deriving asymptotic expansions for large n and m. 



Date: May 1 1, 2010. 

2000 Mathematics Subject Classification. 05A15,60F05. 

Key words and phrases. Urn models, OK Corral, Sampling without replacement, sequential urn model, Jacobi theta function. 
This work was supported by the Austrian Science Foundation FWF, S9608-N13. 

1 



2 



M. KUBA 



1.2. Connection to ordinary Polya-Eggenberger urn models. The dynamics of ordinary Polya-Eggenberger urn 
models are very similar to the dynamics of the two previously descirbed urn models. At the beginning, the urn 
contains n white and m black balls. At every step, we choose a ball at random from the urn - a white ball is drawn 
with probability n/(n + rn) and a black ball is drawn with probability raj (n + m) - and examine its color and put 
it back into the urn and then add/remove balls according to its color by the following rules. If the ball is white, then 
we put a white and b black balls into the urn, while if the ball is black, then c white balls and d black balls are put 
into the urn. The values a, b, c, d € Z are fixed integer values and the urn model is specified by the transition matrix 
AI = i y a c h d )- We can specialise the sequences A = (a„) ne n and B = (/3 m ) m( =N in order to obtain two families of 
Polya-Eggenberger urn models as special cases. 

1.3. Sampling without replacement and generalizations. This classical example corresponds to the urn with tran- 
sition matrix M = ( "q 1 _j ) . In this model, balls are drawn one after another from an um containing balls of two 
different colors and not replaced. What is the probability that k balls of one color remain when balls of the other color 
are all removed? The generalized sampling without replacement urn is specified by ball transition matrix given by 
M = ( with a, d G N. In the general model one asks for the number of white balls when all black balls have 
been removed, starting with a ■ n white balls and d ■ m black balls. Let Y a n.dm denote the random variable counting 
the number of white balls when all black balls have been removed. The random variable Yan.dm is related to X n>m of 
urn model I in the following way. 

Y an ,dm _ y for the special choices A = (a n ) n(i rj = (a • n)„ 6N , B = (p m ) m eN = (d ■ m) meN . 

a 

The most basic case a = d = 1 can be treated by several different methods such as generating functions, lattice path 
counting arguments, etc. The distribution of Y n m is a folklore result and given as follows. 

/n+m— 1 — k\ 
P{Yn,rn = fc} = ' "mN ' , < k < 71. 



1 .4. The OK Corral problem and generalizations. The so-called OK Corral urn serves as a mathematical model of 
the historical gun fight at the OK Corral. This problem was introduced by Williams and Mcllroy in [22] and studied 
recently by several authors using different approaches, leading to very deep and interesting results; see [20, 10, 11, 
12, 19, 3]. Also the urn corresponding to the OK corral problem can be viewed as a basic model in the mathematical 
theory of warfare and conflicts; see [11, 12]. An interpretation is given as follows. Two groups of gunmen, group A 
and group B (with n and m gunmen, respectively), face each other. At every discrete time step, one gunman is chosen 
uniformly at random who then shoots and kills exactly one gunman of the other group. The gunfight ends when one 
group gets completely "eliminated". Several questions are of interest: what is the probability that group A (group B) 
survives, and what is the probability that the gunfight ends with k survivors of group A (group B)? Apparently, the 

described process corresponds to a Polya-Eggenberger urn model with transition matrix M = ( q ) . This model 
was analyzed by Williams and Mcllroy [22], who obtained an interesting result for the expected value of the number 
of survivors. Using martingale arguments and the method of moments Kingman [10] gave limiting distribution results 
for the OK Corral urn model for the total number of survivors. Moreover, Kingman [11] obtained further results in 
a very general setting of Lanchester's theory of warfare. Kingman and Volkov [12] gave a more detailed analysis of 
the balanced OK Corral urn model using an ingenious connection to the famous Friedman urn model; amongst others, 
they derived an explicit result for the number of survivors. In his Ph. D. thesis [19] Puyhaubert, guided by Flajolet, 
extended the results of [10, 12] on the balanced OK Corral um model using analytic combinatoric methods concerning 
the number of survivors of a certain group. His study is based on the connection to the Friedman urn showed in [12]. 
He obtained explicit expression for the probability distribution, the moments, and also reobtained (and refined) most of 
the limiting distribution results reported earlier. Some results of [19] where reported in the work of Flajolet et al. [3], 
see also Flajolet [2]. Stadje [20] obtained several limiting distribution results for the generalized OK Corral urn, as 
introduced below, and also for related urn models with more general transition probabilities. 

A natural generalization of the OK Corral urn model is to consider an um model with transition matrix given by 
M = {° c ~o b ), b,c e N, which may serve as a simple mathematical model of warfare with unequal fire power. 
Let Y cnj b m denote the random variable counting the number of white balls when all black balls have been removed, 
starting with c • n white balls and b ■ m black balls. The random variable Y cn ^ bm is related to the random variable X n>m 
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of urn model II in the following way. 

cn,bm _ x for the special choices A = (a n ) neN = (c • n) neN , B = (/3 m ) m6N = (b ■ m) m£N . 

c 

Concerning the special case c = b = 1 Puyhaubert and Flajolet obtained the following result for the distribution of 

Y 

I 71,771' 

1.5. Notation. Throughout this work we use the notations N = {1, 2, 3, . . . } and No = {0, 1,2,...}. Furthermore, 
we denote with i the imaginary unit, i 2 = —1. We use the abbreviations x- := x(x — 1) • • • (a: — s + 1) for the 
falling factorials. Further we denote with [?] the unsigned Stirling numbers of the first kind and with {^} the Stirling 

(d) 

numbers of the second kind. We use the notation X n — > X for the weak convergence, i. e. the convergence in 
distribution, of the sequence of random variables X n to a random variable X. 

1.6. Overview. We derive explicit formulas for the probability mass functions of the two urn models I and II for 
two given sequences (A, B) of positive numbers A = (a n ) ne ps and B = (/3 m ) m £N assuming that both sequences 
(A, B) satisfy the conditions atj ^ ag and f3j ^ f3g , 1 < j < £ < oo. Our explicit results concerning the distribution 
of X n<m and the resulting expressions for the Polya-Eggenberger urn models generalize the corresponding results of 
[3, 8]. Furthermore, we also study higher dimensional urn models with r > 2 types of balls. Our approach consists 
of constructing a family of functions (/ n ,m(*))n,m£Ni which can be related to the probability generating function of 
the random variables X n>m . It is based on the work of Stadje [20], see also Kingman and Volkov [12]. Stadje did 
derive an integral representation for the probabilites ¥{X n m = k} for urn model II. In this work we revisit Stadje's 
approach and extend it to urn model I in order to obtain integral representations for both urn models, which in turn 
will be used to derive the explicit expressions for the probability mass function P{X n , m = k}. 

Furthermore, concerning the Polya-Eggenberger versions of urn models I and II we will derive various expressions for 
the moments of the random variables X n , m (or more precisely Y an ^m and Y cn ^ m ). Later, we extend urn models I and 

II to higher dimensions by considering r > 2 different types of balls. Moreover, we will prove that both urn models are 
dual to each other in a certain sense, see Section 6, which also extends to the higher dimensional urn models. Finally, 
in Section 7 we study urn model I with weight sequences (a„)„ S N = (n) n( z^ and (/3 m ) m gN = (w 2 )raeM an d show 
that interesting limiting distributions arise, curiously involving the Jacobi Theta function Q(q) = Sriez(~l) n ?™ • 



2. Main results 

2.1. Results for urn model L 

Theorem 1. The probability mass function of the random variable X n _ m , counting the number of remaining white 
balls when all black balls have been drawn in urn model I with weight sequences A — (a„)„ S N, B = (P m )m£N 
(Sampling with replacement with general weights)), is for n, m > 1 and n > k > given by the explicit formula 



m 71 ill 

nX n ,m = k) = ( IJ Ph) ( II a ")E 

h=l h=k+l 1=1 



1 



Jii 1 £=i ( uu (aj + Pii) ( ngi (a -&)) 

(n&)( II a »)Ej=r- ( ,w nm jr- : 

h=i h=k+i e=k yllj=k(aj - at) ) ylU=i(Pi + at) 



assuming that aj 7^ ctg and f3j fig, I < j < £ < 00, and that a a = 0. 



Adapting the weights sequences (A, B) to the special case A = (a ■ n)„ e pf and B = (d ■ m) me N, with a, b E N, 

-a \ 
-dj' 



corresponding to the generalized sampling Polya-Eggenberger urn model with ball transition matrix M = 1 



leads to the following result. 



4 



M. KUBA 



Corollary 1. The probability mass function of the random variable Y an ^. m , counting the number of remaining white 
balls when all black balls have been removed in the generalized sampling without replacement Polya-Eggenberger 
urn model with ball transition matrix M = ( "q" _° d ), starting with a ■ n white and d ■ n black balls, is given by 

P{Y an , dm = ak} = E(-l)'- l U ; B+ «* j = E(- 1 )^ fc M^' 0<k<n. 

1=1 \ n") l=k { m d ) 

The moments of the normalized random variable Y an ^ m — Y an ^ m /a are given by the following expressions. 



/ 3=0 ^ JJ { m d ) 

where {^} denote the Stirling number of the second kind, and x- — x(x — 1) (x — 2) ... (x — (s — 1)) the falling 
factorials. 

2.2. Results for urn model II. 

Theorem 2. The probability mass function of the random variable X n m , counting the number of remaining white 
balls when all black balls have been drawn in urn model II with weight sequences A — (a n ) nl =iq, B = (/3 m ) m gN ( OK 
Corral urn model with general weights), is for n,m> 1 and n > k > 1 given by the explicit formula 



P{X„, m = k} = a k J2 



m+n—k—1 
3 



j=k I Y[l=k{Uj -Otl))[ n h =l("j + Ph) 

m nn-\-m — l — k 
_ \ ~" Pi 

~ ak h ( u;= k Wi + ( n|=i(A - h) 

assuming that aj ^ and j3j ^= /3g, I < j < £ < oo. For n, m > 1 and k — we obtain 

v n+m— 1 

P{Xn,m = 0} = 1 - E 



(nf=i(«i-a.))(nr= 1 K+^)) 



1 



= E 



fe(n?=i(A+«i))(n^j(A-^)) 



Adapting the weights sequences (A, £?) to the special case A = (c • rt) nS N, B = (b ■ m) mS N, with 6, c 6 N corre- 
sponding to the generalized OK-Corral Polya-Eggenberger urn model with ball transition matrix M = (_° c we 
obtain the following results. 

Corollary 2. The probability mass function of the random variable Y C n.bm, counting the number of remaining white 
balls when all black balls have been removed in the generalized OK-Corral Polya-Eggenberger urn model with ball 
transition matrix M = (_° cQ b ) starting with c ■ n white and b ■ n black balls, is given by 

t. h n-k m l™- 1 ) 

ply t - rk\ - W IV"-' W-W m+m-l-k 

(n — K + 1)!(to — lj! C n K f— ' "+;' 

l—l \n— k+l) 



h m " (n—k\ 

' ^(_l)«-'A^r+"- 1 - fc /or l<k<n, 



(n — k)\m\ b 

,„ ... . ... ^, 

P{^c„ bra = 0} = -77 ^— y(-l)" l - / ^4^/" +m - 1 , for fc = 

1 ' ; n\ to- 1)! c n ^ ' ' J 



1=0 ' 



n\(m — IV. - 

K ' i=i . „ - 

Using the explicit expressions for the probabilities and an important result of Puyhaubert [19] and Flajolet [2], origi- 
nally developed for the analyis of the case b = c, we can give a precise description of the moments. 
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Corollary 3. The moments of the random variable Y cn ^ m = Y cn j, m /c, where Y cn ^ m counts the number of remaining 
white balls when all black balls have been removed in the generalized OK-Corral Polya-Eggenberger urn model with 
ball transition matrix M — (® c ~ b ) starting with c ■ n white and b ■ n black balls, are given by the exact formula 

n rn+m\ fm+l\ 

where f n (u) — n\[z n ]F(z, u), g n {u) — n\[z n ]G(z, u), with 

F(z,u) = e u ( e ~ Z+z - 1 ), G(z,u) = e u ( e ~ Z+z - 1 ) f" ue" 1 dt, 

Jo 

and Q(n) = Y^i=o TP denoting Ramanujan's Q-function. Moreover, there exists monic polynomials (M s ) sS n of 
degree 2s such that 



E(M s (Y cn>bm )) 



los « (n+m\ (m+£\ , n fn\ 

/-^i n\m\ b m ^— ' 



K ' £=0 Km) e=o Km) 

where the coefficients s , with 1 < i < 2s, of M s are specified by the triangular system of linear equations 

2s 2s 

^m i , s / i+1 (X) = 0, J2 m ^9i+i(X) = s\2 s X s+1 . 

i=l i=l 

3. Deriving the explicit results: Urn model I 
By definition of urn model I the probabilities ¥{X nm — k} satisfy the following system of recurrence relations. 

a B 
¥{X n . m = k}= ^—¥{X n ^ m = k}+ l m ¥{X n . m -i = fc}, m > 1, n > 1, n > k > 

a n + p m ce n + B m (2) 

¥{X nfi = k} = S k , n , n > 1, k > 0, and P{X , m = fc} = 4,o, m > 1, k > 0. 

In order to prove the result above we first derive a class of formal solution of the system of recurrence relations (2). 
Then, out of the class of formal solutions we determine the proper solution by adapting to the initial conditions 
¥{X n ,o — n} — 1 for n > 0. We will present the proof using a sequence of lemmata. 

Lemma 1. The quantities f n ,m.k{t) G C, defined as 



fn,m,k (^) 



(n^(i-4))(nr =1 (i+^ 



with parameter G C, satisfy for each fixed t G % -K ant/ n, m > 1 and k > 1 a recurrence relation of the 

type (2). Moreover, f n ,m,k{t) W fl meromorphic function in C, vv/f/; /7oZei af f/ze values {—Be | 1 < < to} one/ 
{«j | A; < j < n}. 

Proof. We observe that for arbitrary n,m,k > 1 

fn-l,m.k(t) = fn,m.k(t)(^ V /n,m-l,fe(*) = /n,m.fe(*)fl + "B - Y 

Consequently, 

," fn-l,m,k{t) H fn,m-l,k(t) = ^ n ' m ' k ^ - ( a n — t + B m + t) = f n .m,k(t)i 

a n + p m a n + B m a n + p m \ / 

which proves that f n . mt k(t) satisfies a recurrence relation of the type (2). □ 

Note that we exclude the case k = 0, since we assume that o.q = 0. Since all of the quantities f n , m ,k(t) S C are 
solutions of a recurrence relation of the type (2), we can derive another family of solutions. 
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Lemma 2. The values F n _ nii k — i^—- J, f n ,m,k(t)dt are for n,m,k £ N formal solutions of the system of recur- 
rence relations (2), with parameter ^gC assumed to be a constant independent oft. For n > k > 1 and m > 1 we 
can obtain F n>m ^ by two different contours of integration, 

F n , m ,k = 7^— [ fn,m,k(t)dt = ^— I f n , m ,k(t)dt, with C\ = 71 U 72, C 2 = 7i U 73, 



2iri lr ' ' 2m 

where ~fx{z) = -iz, -R < z < R, <y 2 (z) = Re itp , n/2 < <p < 3vr/2, 7 3 (z) = Re 1 ^-^, n/2 < ip < 3tt/2, 5mc/z 
that R > maxi<e< mt k<j<n{Pti a j}- 

Proof. By definition of f n ,m,k{t) the first stated asseration easily follows. By Lemma 1 f n . m .k(t) is a meromorphic 
function in C, with simple poles at the values {— /3t | 1 < £ < m} and {aj \ k < j < n). By taking the limit R — > oo 
in the two different contours of integration we observe that the contributions of 72 and 73 are zero in the limit, and 
only the contribution of 71 remains. Hence, 

1 r°° ir ir 

fn,m,k(t)dt = - — : / fn,m,k(t)dt = - — ; / fn,m,k{t)dt. 



r\ • I j 1v.1iL.ru \ / . 1 j 1 l . 1 1 L.ru \ / 

2tt« J- l0 o 2tt2 J Cl 2m _ , _ 

□ 

Next we will derive an alternative representations for F n m ^ using the formal residue calculus. 

Lemma 3. Assume that the weight sequence A — (a n ) n< =fn satisfies a, 7^ aj for all i,j € N with i ^ j. Then the 
value F n>mt k admits for n > k > 1 and m, > 1 the following representation. 

m n n . 

F n , m , k = &( n a) ( n «0 e 7^ — .w nm — . v 

h=i h=fe fcfe 1 Ij-fc (aj - ar<) J ^ 1 I i=1 (ft +a,)j 

A.siMme that the weight sequence B = (/3 m ) m gN satisfies Pi ^ /3j for all i, j £ N with i ^ j. Then the value F n ^ m ^ 
admits for n > k > 1 and m, > 1 the following representation. 

1 



111 n £i y ( n;=, + a)) ( n|i (a - a)) 



Proof. By definition of F nm ^ and Lemma 2 we have 



"n.m,k — . / Jn,m,k\t )ClL — 

2m r „ 2m 



n;u(i-i))(nr =1 (i + i) 

By the residue theorem with respect to the poles — aj, k < £ < j, we obtain the representation 

71 m n n ^ 

F n ,m,k = ^ ReS (/™.™. fc ' a ' 1 ) = II Ph) ( \\ a h) ^2 



h=k h=\ h=k e=k ( U7=k(a 3 - a e )j ^ J]™ i(A + 

by the assumption on the weight sequence A. A similar argument leads to the second representation for F n m ]~. □ 

Note that one can alternatively convert between the two expression using the partial fraction identity 

n 

i = y i (3) 

Kj=k( a 3 + x ) ^ k ( x + a h) Uj=k(aj - a h ) 

For example, applying the partial fraction identity to the second expression gives 

ITi Tl Tl 771 

F n , m ,k = ik ( n h) ( n «o e 7 — 1 — v e t 1 

V h=i ' V h=k ' U ( Ul=k (a d - a h )) ti(f3i + a h )( ]jT=i {fit ~ Pi. 
Another application of the partial fraction identity provides the first expression. 
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We still have to adapt the parameter ^ in the representation of F niTn & to the initial conditions P{X„ — k} — 8 nk 
for n,k > 0, and the check the initial conditions. First we will consider the special value F nim>m in order to choose 
the parameter ^ in the representation of F n>rrii k based on the explicit expression. Then we will check if the initial 
conditions are fulfilled. Let k = n. We have 

mm _^ 

V h=i 1 tl (a n + ft) ( n™ 1 (ft - ft) ) 
Application of (3) with respect to the case k = n gives 

11/*=! (ft* +«n) 

Setting m = Owe obtain F n ^, n = £, n a n . In order to adapt to the special cases of the initial conditions P{X„o = 
n} = 1 for n > 1 we choose = l/a^ for n > 1 and k > 1. Then we extend the definition of f n .m,k{t) and F n ^ m ^ 
to cope with the additional boundary cases n = 0, m = and k = 0. We obtain the desired representation by a slight 
modification of the definition of f n ,m,k{t) in the case of k = 0. 

Lemma 4. For n, to > 1 and k > the values P{X n _ m = k} = ^ J c f n ,m,k(t)dt, with 




for k>l, 



t x-7 V' f or fc = ' 

k (-t)(^n?= 1 (i-^)J(reii(i+^)J 

where C 3 = j[ U i 2 U 7^ U 74, 7^(3) = -iz, -i? < z < e, 7^(2) = ee ^ 37r / 2 ^, < 95 < tt, = iz, 

e < z < R and 74(2) = Re l ^^ v \ tt/2 < <p < 3tt/2, such that R > max.i<g< m ^<j< n {l3e, &j} and e < 
{ft, ctj}, are f/ie proper solution of the system of recurrence relations (2) fulfilling the stated initial 
conditions involding the boundary cases n = and m = 0. 



Proof. First we check that the initial conditions P{X n ,o = k} — 5k. n for 77 > 1 and k > 0. For k = n we can use our 
previous computation which lead to the choice = l/ct^, and the stated result easily follows using (3). For k 7^ n 
we use residue theorem and then the partial fraction decomposition 



at 



(x + a k )(x + a k+1 )...(x + a n ) ^ (<r + a*)(n?=*(«i " «*)Y 

Taking the limit .t — > in the identity above provides the needed result. Finally, we check the initial conditions 
P{^o,m — k} — Sk.o for m E N and k <E No- For /c = we immediately obtain ¥{Xo. m = 0} = 1, due to the 
additional simple pole at t = 0. Concerning the remaining cases k > 1 we rely on the integral representation of the 
probabilities J Cs f n . m .k(t)dt. By Lemma 1 f n . m ^ k (t) has no singularities with negative real part (more precisely, 
it has removeable singularities at {aj \ 1 < j < n}). Hence, by the Cauchy integral theorem ¥{Xo, m = k} — 
2^ fc 3 km,k (t)dt = for k > 1 and n = 0. □ 

Remark 1. It seems difficult prove that the initial conditions ¥{Xo, m — k} = Sk.o f° r m E N and fc S No are 
satisfied relying on the other integral representation, or solely on the two explicit expressions valid for n > k > 1 and 

771 > 1. 



Proof of Theorem 1. By Lemma 3 we obtain for m, n > 1 and k > 1 the stated explicit expressions. Our subsequent 
calculations show that stated expressions are also valid for k = 0. By Lemma 1 and Lemma 4 the recurrence relations 
and the initial conditions are are fulfilled. It remains to prove that the stated formulas are also valid in the case k = 0. 
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By the residue theorem we obtain 

(nr^Xn^) 



P{X n>ro - 0} = — / f n ,mfi(t)dt = 1 - V — , ^ ^7 ^ 

27ri ^ fe«*(i%i(«i-«*))(n£i(A+«<) 
nwA)(iL 



E 



(n"=o(«i-^))(n™i(A 
(nr =1 A)(nLx^ 



+ . 



(n|=o(ay -"<))( 112=1 (A + 

subject to a = 0, which proves the first part of the stated result for k — 0. The second formula follows easily using 
the partial fraction decomposition (3). Alternatively, one may change the contour of integration C3 and collect the 
residues at the poles — p\, . . . , — f3 m . □ 

Proof of Corollary 1. In order to obtain the results for the Polya-Eggenberger generalized sampling urn model with 
ball transition matrix M = ( ® d ) , we simple adapt the weights sequences (A, B) to the special case A = (a ■ n)„ e N 
and B = (d ■ m) mg N, with a, b <E N. Concerning the moments we use the well know, see e.g. [7], basic identities 

^ fx + k\ fx + n + l\ ^ fn\ (-l) fc 1 
t„l fc J V n /' fr W^ + fc x(T)' 
Let Y a n,dm denote the scaled random variable Y an ^ m /a. The factorial moments of Y a n,dm stre derived as follows. 

n n m ( m )( k ~ 1+ ^)k- 

nyL, d m) = E fc£p i^„, rfm = ka} = EE(- 1 )'" L+k 

fe=0 t=s <=1 ( tl ° j 

By interchanging summation and using the first identity stated in (4) we get further 

m tm\ (n+s+^\( _ 1 1 Id \s_ ra / m \ 

m n , d j = Ec- 1 )'" 1 1) w -J.L a = ^Et- 1 )'- 1 u 



s + 



The application of the second identity stated in (4) and the fact that 



k=l 

where |?) denotes the Stirling numbers of the second kind lead obtain the following results. 



Km) j=0 > 



\ rn ) 



□ 



4. Deriving the explicit results: Urn model II 

The derivation of the explicit results for urn model II is similar to are previous approach to urn model I; therefore we 
will be more brief. By definition the probabilities P{X„ m = k} satisfy the following system of recurrence relations. 

HXn.ra = &} = „ P{* B -1, TO = k} + P{X n , m _i = A}, m > 1, 71 > 1, U > fc > 

a n + Pra On + Pro (6) 

P{X„, = k} = 4,n, and P{X , m = A;} = 4.o, m > 1. 

As before, we first derive a class of formal solution of the system of recurrence relations (6). Then, out of the class of 
formal solutions we determine the proper solution by adapting to the initial conditions P{X n .0 = n} = 1 for n > 0. 
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Lemma 5. The quantities g nm k(t) G C, defined as 



(n? =1 (i-«j<))(n^i(i+A<))' 

with parameter £fc € C, satisfy for each fixed t G i • M one/ n, m > 1 an<i fc > 1 a recurrence relation of the 
type (6). Moreover, g n ,m,k{t) is a meromorphic function in C, w/f/i poles at the values {— 1/p^ | 1 < I < m} and 
{1/oj \k<j<n}. 

Proof. We observe that for arbitrary n, m, k > 1 

9n-l,m,k(t) — g n ,m,k{t) (l ~ a n tj , g n ,m~ l,k (t) = <?n,m,fc(i) (l + /3 m t) ■ 

Consequently, 



gn-l,m,k{t) H ," a gra.m-l.feft) — ^ n,m ' k ifim ~ Pm,Ol n t + a n + 0> n f3 m t) — g n .m,k(t), 

Ct„ + Pm Ot„. + Dm \ ' 



C^n ~r Pm C^n T Pm ttn T" Ptj 

which proves that g n , m .k(t) satisfies a recurrence relation of the type (2). □ 
Next we derive another family of solutions. 

Lemma 6. The value G n .m,fc = f_°° gn,m,k(t)dt is for n, to, k G N a formal solution of the system of recurrence 
relations (6), w/f/z parameter G C assumed to be a constant independent oft. For n > k > 1 ant/ to > 1 we can 
obtain G n ^ m ,k by two different contours of integration, 

Gn,m,fc = 7T— I 9n,m,k{t)dt = ^— [ ffn,m,fc(*)d*, Ci = 7l U 72, C 2 = 7l U 73, 



27T* 7c 1 ' ' ' 2ni Jc 2 

where ji{z) = -iz, —R < z < R, 72(2) = Re itp , tt/2 < <p < 3vr/2, 73(2) = Re^-v), tt/2 < tp < 3tt/2, such 
that R > maxi<e< m ,k<j<n{l/Pe, 

Proof. By definition of g n ,m,k{t) the first stated asseration easily follows. By Lemma 2 g n ,m,k(t) is a meromorphic 
function in C, with simple poles at the values {— Iffy | 1 < I < to} and {1/ay \ k < j < n}. By taking the limit 
R — s- 00 in the two different contours of integration we observe that the contributions of 72 and 73 are zero in the 
limit, and only the contribution of 71 remains. Hence, 

1 f io ° If If 

^— / g n ,m,k(t)dt = — / g n .m.k{t)dt = — I g n , m ,k(t)dt. 
Z'Ki J _ ioo Zm J Cl l-Ki J C2 

□ 

Next we will derive an alternative representations for G nm ^ using the formal residue calculus. 

Lemma 7. Assume that the weight sequence A — (a n ) ne iq satisfies on 7^ otj for all i,j G N with i 7^ j. Then the 
value G„, m ^ admits for n > k > 1 andm, > 1 the following representation. 

n m+n—k—1 

E j 
7 n \7 m ~ 

3=k [ lVi=k(aj - ai)) ^ UhU( a j + Ph) y 

Assume that the weight sequence B = (J3 m ) m £N satisfies Pi 7^ f3j for all i, j G N with i 7^ j. Then the value G„. m , 
admits for n > k > 1 and to, > 1 the following representation. 

Pi 



Gn,m,k ^ ^ 



fe(n-= k (A+«j))(n^03i-/3fc))" 



Proof. By definition of G nm ^f. and Lemma 2 we obtain the stated results by the residue theorem. □ 

We adapt the parameter £k in the representation of G n , m ,k to the initial conditions P{X„.o = k} = 5„.k for n > 1 
and k > and obtain after a few calculations the necessary condition £ fe = for k > 1. 
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Lemma 8. For n, m > 1 and fc > f/ze values ¥{X n m = fc} = ^ g n ,m,k(f)dt, with 



w V' 



i 



(-*) ( n?=i(i-a^)) ( n£i(i+&*)) 



/or fc = o, 



w/zer<? C 3 = 7^ U 7^ U 7^ U 74, 7^) = -iz, -R < z < e, j' 2 (z) = £e ^/^-v) > < ip < n, ^(z) = iz, 
e < z < R and j' 4 (z) — Re i ^ T ~ lf '\ n/2 < tp < 3n/2, such that R > maxi<e<m,k<j<n{^//3e, l/ a j} an d £ < 
minK^< m ,k<j< n {l/ Pi, l/aj}, are the proper solution of the system of recurrence relations (2) fulfilling the stated 
initial conditions involding the boundary cases n = and m = 0. 

Proof. The initial conditions PjX^o = k} = 5k, n f° r n > 1 and fc > 1 are satisfied since = a^. Finally, we check 
the initial conditions ¥{Xo >m = fc} = 6k,o for to e N and fc 6 No. For fc = we otain P{Xo, TO = 0} = 1 due to the 
additional simple pole at t = 0. Concerning the remaining cases fc > 1 we rely on the integral representation of the 
probabilities ^ f c g n ,m,k{t)dt. By Lemma 1 g n ,m,k(t) has no singularities with negative real part (more precisely, 
it has removeable singularities at {1/oy | 1 < j < n}). Hence, by the Cauchy integral theorem P{Xo, m = fc} = 
2¥i Ic 3 fa,m,k(t)dt = for fc > 1 and n = 0. □ 

Finally, the probabilities F{X n m — 0}, with n, m > are calculated in two different ways using the residue theorem. 

Proof of Corollary 2. In order to obtain the results for the Polya-Eggenberger generalized sampling urn model with 
ball transition matrix M = ( _° c ~ b ), we adapt the weights sequences (A, B) to the special case A — (c ■ n)„ S N and 
B = (b ■ TO,) me jj, with a, b € N. Concerning the moments of Y cn ,bm — Ycnbm/ C we proceed as follows. We obtain 
the ordinary moments E(l^ bm ) with s > 1, using the explicit expressions for the probabilities P{Y cn ^ m = cfc} by 
summation. 

n 71 1 s+1 m n in—k\ 

E(y c ;, bm ) = E ^ = cfc > = E T^yw ^ E(-D^ 7 M^ + "- 1 - fc - 

fe=l fc=l v ; ' fcO I m I 

Interchanging the order of summation leads to 



C"° x pn+n £ ^ k s ^t~ 



E { Y cn,bm) - , m E^ 1 )" £ / m+ ^N / ^^"X] 



fc-1 



v_ 1 «- t « = fM^ r +ny fc!r fe fc s . 

The second sum was studied in detail in [19]. We will use the following important result. 
Lemma 9 (Puyhaubert [19]). Let the functions F(z, u) and W(z, u) be defined as 

F{z, u) = e w ( e "" + *- 1 ) , G(z, u) = e"( e " 2+ ^ 1 ) [' ue-'e^i 6 "^ 1 ) dt. 



Jo 

Furthermore, let f n (u) = n\[z n ]F(z, u) and g n (u) = nl[z n ]G(z, u), where the polynomial f n has degree [f J and 
g n has degree L^^J- ^ ne ^ ea ding coefficient of f2 n is given by (2n — 1)!! = (2n — l)(2n — 3) • ■ ■ 3 and the leading 
coefficient o/^n+i is ^iven fey (2n)!! = (2n)(2n — 2) ■ • • 2. Moreover, let Q denote Ramanujan's Q-function defined 
as 

n n(n — 1) n\ ^— * n 1 



Q(n = 1 + - + V 2 7 +••• + — = E 

i=0 

The sum 2~2k=i (kZ.\)k^ k k s can be represented in terms of f n , g n and Q as follows. 

E ( e k Z\) k - rkkS = \{fs + x{t)Q{i)+9s + im- 

k=l ^ ' 
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The result above is proven in Puyhaubert [19], p. 117-118, we refer the reader to the proof therein. As pointed out 
in [19] the function F(z,u) is closely related to Stirling numbers of the second kind. Note that in [19] there is a 
tiny misprint, namely the factor n\ in the definition of /„ and g n is missing. For a comprehensive discussion of 
Ramanujan's Q-function we refer the reader to the work of Flajolet et al.[5]. Using Lemma 9 we immediately obtain 
the stated description of the moments of of Y C n,bm- Furthermore, the second assertion follows easily by comparison 
of the degrees in f n (X) and g n (X), similar to [19]. The second expression for E(M s (Y crit t, m )) is easily obtained by 
simple manipulations of the binomial coefficients. □ 



5. Extension to higher dimensional models 



It is natural to ask whether the simple explicit expressions for the probabilities ¥{X n ^ m = k}, as given in Theorem 1, 
can be extended to higher dimensional models, i.e. involving more than two balls. In the following we will generalize 
urn models I and II by considering urns with r > 2 different types of balls. 

At the beginning, the urn contains ri£ balls of type £ balls, with ne > 1 for 1 < t < r. Subsequently, we will 
use the multiindex notation n = (n 1: . . . ,n r ). Associated to the urn model are r sequences of positive numbers 
A$ — (c4i]) n< eN> for 1 < £ < r. At every step, we draw a ball from the urn according to the numbers of balls 
present in the urn with respect to the sequences {A^)i<£< r , subject to the two models defined below. The choosen 
ball is discarded and the sampling procedure continues until type r of balls is completely drawn. We are interested in 
the joint distribution of type 1 up to type r — 1 balls when all type r balls have been drawn.This corresponds again to 
the absorbing region A = A\ U A%, with A\ = {(«i, . . . , n r _i, 0) : rii, . . . , rv_i > 0} and _4 2 = {(0, ■ ■ • , 0, n r ) : 
n r > 1} We study the random vector X n = (X^\ . . . , X„ where = Xn\,...,n r counts the number of type £ 
balls, starting with n balls of types 1, . . . , r, i. e. ri£ balls of £ units, 1 < £ < r. Subsequently, we will derive explicit 
expressions for the probabilities P{X n = k}, for k = (k\ , . . . , fc r _i), with kj > 0, for the following two urn models, 
extending our previous results corresponding to the case r = 2. 



• Urn model I. Assume that the urn contains ri( balls of type £ balls, 1 < £ < r, with n r > and at least one of 
the rij greater than zero, 1 < j < r— 1. A ball of type £ is drawn with probability ctn]/ Y^j=i 1 < £ < ~r, 

with the additional assumption that a = for 1 < £ < r. 

• Urn model II. Assume that the urn contains ng balls of type £ balls, 1 < i < r, with n r > and at 
least one of the rij greater than zero, 1 < j < r — 1. A ball of type £ is drawn with probability (1 — 

^,o)m=i(^) 1 ^'7ELi(i-^ h ,o)m=i(^) 1 ^'°- 

Note that urn model I is the natural extension of the urn model I with two types balls. By suitable choice of the 

sequences A$ = (ajf])„^gN = ( a e ■ i<)n,£N, for 1 < £ < r, we obtain the standard Polya-Eggenberger urn with r 
different kind of balls, and ball transition matrix M given by 



M = 



V 



\ 



7 



(7) 



In contrast the extension of urn model II is non-standard in the sense that for any specification of the sequences 
= (c4f])n £ GN, 1 < I < t, we do not obtain ordinary higher dimensional Polya-Eggenberger urn models anymore. 
The reason for this will become clear in Section 6, where we uncover the duality relation between these two urn 
models. 



5.1. Explicit results for the probabilities. 

Theorem 3. The distribution of the random vector X n , counting the number of remaining type 1, . . . , r — 1 balls when 
all type r balls have been drawn in urn model I, is for rij > 1 and rij > kj > 0, with 1 < j < r, given by the explicit 
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formula 

;n^i«? ] )nS(n5^ + i«e) 



[X n = k}=£-. E 



1 — fcl fc 



assuming that off ^ a 1 ^ , 1 < h < £ < oo, 1 < j < r. 



Corollary 4. The distribution of the random vector Y an = (Yan > ■ ■ - j ^in counting the number of type £ balls, 
1 < £ < r — 1, in the sampling without replacement Polya-Eggenberger urn model with r different types of balls (J), 
when all type r balls have been drawn, starting with CLjfij balls of type j, 1 < j < r, is given by 

rr-l 



(Y an =ak } = y: ■■■ e - ll • 



i=fci £,._i=fe r ._i 



The mixed factorial moments El YHj=i (if 11 '"'') ) of Y an are given by the exact formula 



E 



Theorem 4. 77ie distribution of the random vector X n , counting the number of remaining type 1, . . . , r — 1 £>a//s w/zen 
a/Z fy/?e r £>aZZ.v /lave been drawn in urn model II, is for rij > 1 and rij > kj > 1, w/f/i 1 < j < r, given by the explicit 
formula 



p{x n = k }= £ ... J2 



n /= i ( n,=i < J + «v e 9= i 3?i J n J=1 n^= fcJ (< - < J ) 



fi=fci f»._i=fc r _: 



assuming that o$ ^ 0$ , 1 < h < £ < oo, 1 < j < r. 



Similar but slightly more involved expressions exist for urn model II in the cases kj = 0, for 1 < j < i — 1. 

5.2. Deriving the explicit results. Subsequently, we will briefly discuss the proof of Theorem 3 and comment on 
the proofs of Corollary 4 of Theorem 4, which are extensions of the proofs of the two dimensional cases, r — 2. By 
definition of urn model I the probabilities P{X n = k} satisfy the following system of recurrence relations for n r > 1, 
nh > 0, 1 < h < r — 1, with at least one of the rih > 0. 

P{X n = k} = nt b1 P{X B - e , = k}. (8) 

e=i l^j e =i a ' n i 

where denotes the £-th unit vector, with respect to the additional assumption that off = for 1 < £ < r. The initial 
conditions are given by 

r-l 

P{^m,...,n r _i,o = k} = JJ 5fc i>nj , rij > 0, kj > 0, with at least one rij > 0, 

(9) 

r— 1 

P{X ,...,0,„ r = k} = n **i,0. «r > 1, % > 0. 

In order to prove the result above we first derive a class of formal solution of the system of recurrence relations. Then, 
out of the class of formal solutions we determine the proper solution by adapting to the initial conditions. The crucial 
part of the proof is the following result. 
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Lemma 10. The quantities / n ,k(t) € C, defined as 

/n,k(t) = 



w/f/j parameter £k(t), satisfy for each fixed tt G i • M anii > 1 ant/ fc^ > 1, /or 1 < I < r — 1, awe/ n r > 1, a 
recurrence relation of the type (8). 

In order to validate Lemma 10 we observe that 

r [I] r-l [£] , [r] v^' - 1 + 

= -Vw^E-eO - +«e0 + ^)) - 

— 1 v £=1 U-n^ t^n r / 

As in the case r = 2 we derive another family of solutions using the Cauchy integral theorem. 

Lemma 11. The values F n ^ = ( 2 iri) r — 1 Ici ' ■ • /c 1 fn,k{tydt r • ■ • dii are for n,m,k £ N formal solutions of the 
system of recurrence relations (8), with parameter £ic assumed to be a constant independent oft, where the interior of 
the contours of integration Ci are containing the the values a£ , 1 < t < r — 1 and kg < hf < ri£, plus the origin, 
and the orientation of these closed curves is mathematically negative (clockwise). 

It still remains to adapt to the initial conditions (9). It turns out that the choice 

&(t) ' 



provides the right solution. Note that £k(t) depends to t, but only in the case of k( — 0, 1 < £ < r — 1. We skip the 
long and involved calculations, which are in principle similar to calculations for r = 2. 

Concerning the result for the Polya-Eggenberger urn model with ball transition matrix M given by (7) we adapt the 
sequences = (an\) ne eN = (at. 1 ni)n t en, for 1 < I < r, and simplify the result of Theorem 3 to obtain the first 

statement of Corollary 4. The mixed factorial moments Ef YYj=i ) °f Y an are derived as follows. 



E 



l — 1 ,, \j] \ m n r ~l r— 1 

n(^ - E nv— *}n«f 

2=1 3 7 ki=0 fc r _i=0 2=1 



f it 1 f™ 3 ) 

e ••■ e l p;; nEH) v , 

E r — ) i=lfe . =fl . V*2 



"' 5 f n ' V) r-l ^ ^ , /£ \ /jfe- 



£l=Sl £ r _ 1= s r _i 



^ v . (in ! o^0w, _ n; 



where 6g tS denotes the Kronecker-delta function. The results for urn model II are derived analogously to the proof of 
Theorem 3 by first deriving a class of formal solution of the system of recurrence relations and then adaptin to the 
initial conditions. 



6. Duality of the two urn models 

One can notice similarities between the explicit results for the probabilities ¥{X n _ m = k} for urn models I and II 
in the case r = 2 and also the similarities of the proofs. However, the higher dimensional models seem to be very 
different. We will show that urn models I and II are closely connected to each other, as they satisfy a certain duality 
relation, and are essentially the same urn model, also for the higher dimensional cases. 
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Theorem 5. Let F{X n>m aa,b,i] — k} denote the probability that k white balls remain when all black balls have 
been drawn in urn model I with weight sequences A — (a n )neN» B — (/3 m )meN an d ^°{X n m rj g = k} the 
corresponding probability in urn model II with weight sequences A — (a n )n£N> B — (/3 m )meN- The probabilities 
F{X ntm! [A,B,l] = k} and V{X n m rj ^ 7J j = k} are dual to each other, i.e. they are related in the following way. 

P{4,m,[A,Bj] = k} = F{X n>m ,a,b,ii] = k }' f or a n = ^~ > Pm = y- , far all n,m £ N. 

a " Pm 

Remark 2. According to the result above, there is only one urn model lurking behind urn models I and II, which can 
be formulated either as an urn model of OK Corral type or of sampling without replacement type with general weight 
sequences. 

Proof. We observe that the recurrence relations, as stated in (2), (6) can be transformed in the following way. 

W{Xn,m,[A,B,I] = k} = W{X n - 1>mi [ AiB ,I] = k } H " P{^n,m-1,[A,.B,J] = k} 

C^n ~r Pm T" Pm 

1 1 

= ~ — ™ i ^{^n-l,m,[A,B,I] = k} H 1 1 ^ {-^n,m-l,[A,S,J] = ^} 

-P{X n _i iTO) [A,B,Z] = fc} + - F{X n<m -i t [ A< B,T\ = k} 



Pm ~r C^n Pm T" Oi n 

^V^n-l,m,[A,B,II\ = k} + ~ _ P{^"n,m-1,[A,B,II] = ^1 



Pm i OL n Pra \ ™n 

= n^.mJA.B,//] = 

Since the initial conditions coincide we can transform one urn model into the other, which proves the stated result. □ 

In particular, the duality relation explains why urn model II looks so different in higher dimensions . 

Theorem 6. Let P{X n ja,j] = k} denote the probability that kj type j balls remain when type r balls have been 

drawn in urn model I with weight sequences A = (A\, . . . , A r ) with Aj = (an] ) nj eN, and P{X n ^ 7J j = k} the 

corresponding probability in urn model II with weight sequences A = (Ax, . . . , A r ) with Aj = (ctn|)n,-gN- The 
probabilities P{X n jA,i] = k} and P{X n ^ ri , = k} are dual to each other, i.e. they are related in the following 
way. 

P{X„, [A ,j] = k} = P{X n [A m = k}, for = —nr, far all n, > 1, with 1 < j < r. 

Proof. We observe that the recurrence relations, as stated in (2), (6) can be transformed in the following way. 

P{X n>[Ai/] = k} = ]T ni tf1 P{X n _ eti[A , J] = k} = £ u "= ia "* P{X n _ ef , [A , f] = k} 



=1 l^je=l a ™j fclT r , ^-n- 



it i a [j] IT i a H 

= v^r Vr T[hT P { X n-e £ ,[AJ] = k } = ^ _r T[S[ P { X n-e^[A,//] = k } 

1=1 2^, = 1 1 lh=l Otn h t =i Z^j=l 1 l£=l 

= F { X „,[A,//] = k}, 

for rij > 0, 1 < j < r. One still has to take into accounts the various cases of rij = 0, 1 < j < r — 1 which leads to 
the Kronecker delta terms in description of urn model two. Since the initial conditions coincide we can transform one 
urn model into the other, which proves the stated result. □ 
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7. Discussion of a special case 

In this section we study in more detail a special case of urn model I with respect to the choices A = (a n ) n eN = 
(n) nS N and B = (f3 m ) me ^ = (m 2 ) me ^ for the weight sequences (A, B). We are interested in the various limiting 
distributions arising for max{n, m} — » oo. Interestingly, the Jacobi Theta function 0(g) = J2nez(~ 1)"<7™ naturally 
arises. 

Theorem 7. The random variable X n _ m , counting the number of remaining white balls when all black balls have 
been drawn in urn model I with weight sequences A = (n) ne ^, B = (m 2 ) me ^, satisfies the following limiting laws. 

(1) Case arbitrary but fixed m G N and n — > oo: X n ^ m /n converges in distribution to a continuous random 
variable Y m with support [0, 1]: 

^ Y m , ¥{Y m <x}= f f m (q)dq, < x < 1, 

with density function f m (q) being given by 

m (m\ 

W9)=2E(-i) M A^ _1 - 

1=1 \ m I 

We have also have convergence of all positive integer moments of X nm , 



:(^-)%E(^), E(Y r s r , 



n ) y mn K m ' r(m + 1 - iy/a)T(m + 1 + i^fs) sMiTT^/i' 

where i denotes the imaginary unit i 2 = —1. 
(2) Case arbitrary but fixed n <G N and m — »■ oo: X n , m converges in distribution to a discrete random variable 

x n , m ^ Zn , F { z„=, } =x:(-i)^f;)ff)^=E(-^- i ^? ! 0<*<n. 

£=k 



We have also have convergence of all positive integer moments of X n ^ r 



(-1) 



3-1. 



sinriTry^J ' 



E(X^ m ) h. E(Z«), E(Z' n ) = EX E { • j 

£=i j=e 

(3) Case min{n, to} — > oo: X n _ m /n converges in distribution to a random variable W with support [0, 1], whose 
distribution function can be expressed in terms of the Jacobi Theta function 9(g) = J2nez(~ 1)™?" . 

^ W; P{l^ < g} = 1 - 6(g), < q < 1. 
n 

We /zave a/so /zave convergence of all positive integer moments of X nm , 



E 



(^)%E(I^ S ), E(W S ) 



sinh 7r-^/s ' 



We can slightly extend the previous theorem by giving a more detailed description of the limit laws by studying the 
limiting distributions of the random variables Y m and Z n . 

Corollary 5. The limit laws of the random variable X n<m , counting the number of remaining white balls when all 
black balls have been drawn in urn model I with weight sequences A = (n) neN , B = (m 2 ) meN , can be described by 
the following diagram concerning the convergence in distribution of the various random variables. 

v n->oo v 

normalization — 



m — > oo 



min{fi,m}->oo 



Z n > w 

normalization ^ 
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Corollary 6. The random variables Y m , with meN, and W of Theorem 7 satify the following decompositions. 



Y m = exp ( ~ E 



£2 J ' 



oo 



where the ei are independent and identically distributed random variables with exponential distribution Exp(l). 

Proofs of Theorem 7 and Corollary 5. In the following we will outline the main steps of the proofs of Theorem 7 and 
Corollary 5. Our starting points are the explicit expressions for P{X nm = k} as given in Theorem 1, 



nXn,m = k} = (™0 2 ^,E TF 



1 



^ixu(^ 2 )ngi(* 2 



i i 



0' 



By routine manipulations we derive the alternative represenation 



¥{X n 



k} 



rn Ik- 



k+e — 1\ (m\ 
k 



ln+£ 2 \ lrn+£\ 
V n J V m / 



e=k 



n 



Concerning case (1) such that m € N is fixed and n — » oo we proceed by deriving a local limit law, i.e. we study 
n ■ V{X n ^ m — k} for k = [n- q\ with q G [0, 1]. After manipulation of the stated products we apply Stirling's formula 
for the Gamma function 



r(x) 

to the quotient of Gamma functions 



/2tt 



(10) 



/fe+r-i\ 

( k I _ j/2 



n n + i)T(k + e) e<f-\ i , 

T{n + l + P)T{k + l) n 1 V j ' 



Thus we obtain a local limit theorem, which in turn implies the convergence in distribution after a few routine manipu- 
lations. In order to prove the moment convergence we study first the factorial moments E(X^ TO ) = M(X n}Tn (X n , m — 
1) . . . (X njn — s + 1)). We obtain a closed form expression in the following way. 



fra\ n 



= ^F{x ntm = k} = 2£(-i)^^5>a-^ 

£=1 \ m I k=0 V n ) 



k+e z -i\ 

k )_ 

2\ ' 



k=0 



with s > 1. We use the summation identities 

X + s + 1 



fe=i 



e p. z a n 



e=i 



1=1 



P 



n 



i 

S + P' 



for s e N, where the second identity easily follows from the partial fraction decomposition (3). We obtain the simple 
result 

m l m \ „ 2 m „ 2 

1=1 Km) 1=1 

Consequently, we can write the factorial moments as polynomials in n of degree s, 

e 



E(X* 



n p2 1 „ e 



i=i i=i 

where the [?] denote the unsigned Stirling numbers of the first kind, also known as Stirling cycle numbers. In order 
to obtain the ordinary moments we proceed as follows. 



x =Erih(^o=Er^n / Av^ \ (-^-^e^e 

h=l J 1=1 L J 1=1 j=l 



3 = 1 



h=l 



h 2 +j' 
(ID 



This implies that for n —> oo the ordinary moments satisfy the expansion 



P 



P 



0(n s 
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which in turn implies the convergence of the normalized moments. Note that 

t\J^ = H) 2 

j_l£ 2 +s T(m + 1 - iy/s)T(m + 1 + iy/s) sinh Try's' 
where i denotes the imaginary unit i 2 = — 1. 

Concerning case (2) such that n € N is fixed and m — > oo we simply observe that the probability mass function 

F{x n , m = k } = ±(-iy- k ( n \ (f) n nz n = k}, 

e=k V / V / J=1 J 

due to 

^2 



f 2 + s m-yoo J- -I- £ 2 + s sinh Tr\/i" 

The moment convergence follows immediately from the explicit expression of the ordinary moments of X n>m . For 
the final case (3) min{m, n} — > oo we check the convergence of the moments using the explicit expression (11) 

'X n . m \ s kJ~s 



E 



> E(W S 

V 71 / mini 7i, m|— >oo 



min{n,m}=foo Smh IX \fs 

By Carleman's criterion the moments define a unique distribution. In order to identitfy the distribution function we 
observe, using Stirling's formula (10), that the probability mass function ¥{X n<m — k} satisfies a local limit law of 
the form 

oo 

nP{X n , m = k} = (2 ^(-ly'-W 3 - 1 ) (1 + o(l)) = &(q)(l + o(l)). 

for k = [n ■ q\ with q e (0, 1), where Q(q) = E„ e z(~ 1 )"9™ 2 = 1 + 2 E„>i(- 1 )"9™ 2 denotes the Jacobi Theta 
function. This implies after a few routine calculations, consisting of deriving the asymptotics of the distribution 
function P{X n m < qn} = Efco ^{-^n,m = the convergence in distribution of the random variable J n m /fi 
to a random variable W with distribution function 1 — <d(q) and support [0, 1]. It is not immediately obvious that 
lim^x- P{W < q} = 1. However, by using the Jacobi Triple-product-identity 

00 

e( q ) = ]T(-i)V l2 = 1 + 2 ]T(-i) V 2 = - 9 2j )(i - q 21 - 1 ) 2 , 

nGZ n>l j = l 

we see that P{W < q} = 1 - 9(g) = 1 - H°° =1 (l - ^X 1 - Q 2 ^ 1 ) 2 , such that 1 - O(q) is indeed a well defined 
distribution function with support [0, 1]. Concerning Corollary 5 we simply observe that 

— — ^— > W, and also F m d > W, 
which can be easily checked from the explicit expressions of the moments. □ 

Proof of Corollary 6. Let e = exp(A). The Laplace transform E(e~ se ) of e is given by E(e~ se ) = j^pr. L et ( e i)ieN 
be i.i.d. Exp(l) distributed random variables. Using the fact that c • e = Exp(i), we obtain 



E 



m m rn 1 

(exp ( - s = l[E(e-^) = J] 7—^ • 

£=1 £=1 i=\ 



Since the moments random variable Y. m and are given by 



£=1 v £=1 

we obtain the stated decompositions, 

m 00 

Y m =exp(-^^) ! W ( = ) exp(-^g 



=1 £=1 

□ 



18 



M. KUBA 



Note that such decompositions are closely related to Brownian excursions and Zeta functions, see the work of Biane, 
Pitman and Yor [1]. 

There exist fully analogous results for closely related weight sequences, also involving functions related to the Jacobi 
Theta function. We want to point out two particular choices, namely A = (a„)neN — ( n )neN> B = {0m)meN = 
(( m ^ ))meNi an d A = (n)n£N, B = ((m — |) 2 ) m eN For example, choosing the beforehand mentioned weight 
sequences A — (n) n£ N and B — (("^jjmeB, the random variable X n , m /n converges for imn{rn,n} —> oo in 
distribution and with convergence of all integer moments to a random variable W with support [0,1], such that 

F{W < q} = 1 -^{-ifVlt + lV^*), with moments E(W S ) = -. ^ , . 

7^o cosh (f V8.s - 1) 

Using an identity of Jacobi for the cube of the Euler function <f>(q) = Yl^Lii^- ~~ <?")> namely 

0^ =£(-1)^+1)/^), 

e>o 

one can observe that ¥{W < q} = 1 — 4> 3 (q) is indeed a well defined distribution function satisfying ¥{W < 1} = 1. 
For the weight sequences A = (n) ne ^ and B — ((m — \) 2 ) m en one encounters a similar random variable W with 
support [0, 1], 

4^ .. 1 



F{W <q} = - V(-l) f ^ , with moments E{W S ) - 

7T ^— ' 21—1 COSh ( 7T 



£>1 



Remark 3. Since the product Ilfci^ r /(^ r + s ) converges for arbitrary real r > 1, we expect that the results of 
Theorem 7 and Corollary 5 can be extended at least to weight sequences of the form (n)„ S N, (m r ) me N, for r > 
1. Moreover, decompositions into sums of exponentially distributed random variables, similar to Corollary 6, are 
expected to appear, which are again closely related to the distributions discussed by Biane, Pitman and Yor [1]. 

8. Conclusion and Outlook 

We have studied two urn models with general weights, generalizing two well known Polya-Eggenberger urn models, 
namely the sampling without replacement urn and the OK Corral urn model. Assuming that the urn contains two types 
of balls, we derived explicit results for the distribution of the number of white balls, when all black have been drawn. 
As a byproduct, we obtained results for the aforehand mentioned Polya-Eggenberger urn models. We also studied 
higher dimensional urn models and managed to obtain again explicit expressions. Furthermore, we have shown that 
the two urn models studied in this paper are essentially the same, since they obey a duality relation, also valid for the 
higher dimensional urn models. Examplarily, we discussed in the present work a particular example of urn model I 
with weight sequences A = (n) n grj, B — (m 2 ) m£ N and obtained a detailed description of the limit laws, curiously 
involving the the Jacobi Theta function <d(q) = J2 n ez(~^) n 1 n • Since both urn models can be handled together, the 
author is currently studying the various limit laws arising in the Sampling without replacement/OK Corral urn models, 
trying to unify the previous results in literature. In particular, the weight sequences A = (n°), ie n, B = (m b ) m6 N, 
with a, b e K \ {0} lead to a variety of different limiting distributions, some of which are closely related to Brownian 
excursions and the exponential constructions of Biane, Pitman and Yor [1]. 
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